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We address theoretically the optical parametric oscillator based on semiconductor cavity exciton- 
polaritons under a pulsed excitation. A "hyperspin" formalism is developed which allows, in the case 
of large number of polaritons, to reduce quantum dynamics of the parametric oscillator wavefunction 
to the Liouville equation for the classical probability distribution. Implications for the statistics of 
polariton ensembles are analyzed. 



I. INTRODUCTION 

Semiconductor microcavities with embedded quantum 
wells exhibit a rich variety of unusual light-matter cou- 
pling effects. In the strong coupling regime excitons 
and photons are mixed into a new kind of quasiparti- 
cles, known as cavity exciton-polaritons 1 . These parti- 
cles inherit sharp energy dispersion of cavity photons and 
strong interaction nonlinearities of excitons. The bosonic 
nature of cavity polaritons allows to observe a num- 
ber of coherent phenomena such as stimulated scatter- 
ing of exciton-polaritons and even, perhaps, their Bose- 
condensationS^. 

A non-parabolic shape of the lower dispersion branch 
of exciton-polaritons having a sharp mininum and flat 
wings allows for the so called parametric process when 
a pair of polaritons from the pump scatter into the sig- 
nal and idler states with both energy and momentum 
conservedM^ (see Figure 1). The parametric scatter- 
ing process was studied in detail experimentally under 
continuous- wave (cw) pu mping in R efs. 07 8 9 and under 
pulsed excitation in Refs. Il0lllll2l Previous theoretical 
treatment was based on the Heisenberg equationsi&A^i^. 
The quantum state of the "pumped" polaritons has been 
described by a classical non-fluctuating field which is 
valid only in the case of cw pumping. This approach al- 
lowed to obtain a closed set of equations, describing para- 
metric amplification. However, effects of the parametric 
oscillations, when due to scattering polaritons from sig- 
nal and idler states return to the pump state are lost in 
this description. Moreover, solutions obtained in Refs. 

become unstable above the stimulation thresh- 
old. The semi-classical treatment based on the approach 
analogous to the methods used for the description of the 
four-wave mixing phenomena presented in Refs. Il6ll7l 
gives no insight into the microscopic effects. The general 
theory of the parametric oscillator, relaxing the above 
approximations, is lacking in the literature, to the best 
of our knowledge. 

The aim of this paper is to develop a general quan- 
tum formalism allowing an analytical treatment of the 
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parametric oscillator under the pulsed pumping, both 
below and above the stimulation threshold, provided the 
total number of polaritons is large. We introduce the 
hyperspin pseudovector whose components describe the 
pair correlations between pump, signal and idler states. 
Starting from the Heisenberg equations for hyperspin 
components and treating them as non-correlated quanti- 
ties we obtain classically stable trajectories for the hype- 
spin dynamics. Then, within an approach close to that 
used earlier by one of us to describe the dynamics of 
large total spins of magnetic polarongi^, we analyze the 
Schroedinger equation in the hyperspin space, and show 
that the number of relevant variables in the equation can 
be reduced, and its quasi-classical solution can be found. 
It turns out that the dynamics of the parametric oscilla- 
tor can be described by the Liouville equation where the 
squared wavefunction of the parametric oscillator plays 
the role of a classical distribution function. We show that 
polaritons pass from the pump to the signal state with 
some delay, and that the populations of the pump, signal 
and idler states demonstrate dumped oscillations about 
their average values. We also obtain values of the second- 
order coherence in the steady-state regime and demon- 
strate that parametric oscillations of polaritons between 
signal, pump and idler lead to disappearance of the initial 
coherence. 

II. THEORY 

We consider a typical experimental situation when the 
polaritons are excited in the lower branch of dispersion 
under a "magic angle" (see Fig. 1). The process of the 
parametric scattring of two pump polaritons into the sig- 
nal and idler states can be described by the following 
Hamiltonian 

H = Hq +Hint, (1) 

where the first term Ho = E s ala s + + E % a\ai 

is the free propagation term, Ek (k = s,p,i) are the 
energies of signal, pump and idler states, and a k and ak 
are bosonic creation and annihilation operators for each 
state. The second term in (Q, describing the polariton- 
polariton interaction, readsi*iiiiii^: 

Hint = V(a p a p a\a\ + a s aia p al), (2) 
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where V is the constant of polariton-polariton interac- 
tion, the first term in brackets describes scattering from 
the pump into the signal and idler, and the last one ac- 
counts for the reverse process. 
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Therefore, we introduce a pseudospin vector I with 
components X, Y, Z and the value / = N/2: 

N 2 N 

X 2 + Y 2 + Z 2 = I{I + l) = — + -, (6) 

where N is the total number of polaritons. 

For further treatment, it is more convenient to use an- 
other set of variables £,r),Z where £ and 77 are defined as 



X + Y 



X-Y 



(7) 







The commutation relations for 77 and £ are the same as 
for X and Y respectively because 77 and £ are obtained as 
a result of rotation of the coordinate frame in the (X, Y) 
plane. In terms of new variables, the interaction part of 
the Hamiltonian Eq. (J5J can be written as 



FIG. 1: Polariton energy (w±) vs incident angle at detuning 
A = (solid). Cavity and exciton energies (ui cav , are 
shown by dashed lines. Two pump polaritons (•) scatter into 
signal (o) and idler (<£>). (Adapted from Ref. y). 



Hint = 2V(^ + 7$ 



(8) 



Using commutation relations (Eqs. (JSJ and Q) we ob- 
tain the following Heisenberg equations for £,77 and Z 
(dot means time derivative) 



A. Degenerate parametric oscillator 



Z= - — {71 



a, 



(9) 



To start with, we consider a simplified 2-level model, 
where signal and idler states coincide [this model was 
invoked earlier to describe experiments on degenerated 4- 
wave mixinglSi] . The interaction part of the Hamiltonian 
for this system can be represented as 



Hint = V(a p a p a\a\ + a s a s a p al). 



(3) 



It is well known that quantum mechanical description of 
two selected states of one particle is possible in the terms 
of a fictitious spin 1/2 (Ref. 1201) . This description can be 
easily expanded to a system of N bosons occupying two 
quantum-mechanical states. Let us introduce operators 
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X = - (a s al + Opal) 



(4) 



2V 

f) = —(Zri + r)Z), 



2V 



Since the number of polaritons is very large, / 3> 1 and 
we can treat the pseudospin as a classical vector and 
consider trajectories of the end of this vector on a sphere 
with the radius / (see Fig. 2). Two meridional circles 
(77 = and £ = 0) are the stable trajectories, i.e. any 
trajectory starting near one of the poles (77 = 0, £ = 
0, Z = I and 7/ = 0, £ = 0, Z = —I) very soon comes very 
close to one of these circles. The equations of motion on 
the meridional circles are exactly solvable. For example, 
for £ = we make a substitution Z = I cos ip and 77 = 
/ sin (p which gives a single equation for ip 



Y = — - (a s a. 



a v al 



AVI 



■ sin ip 



(10) 



Z =0 i a l a P - a l a s) 



whose mean values give the difference of occupation num- 
bers of the two states (Z) and second-order correlations 
between them {X and Y) . These operators obey the com- 
mutation relations of angular momentum components: 



[X,Y]=iZ, [Z,X]=iY, %Z) = iX 



(5) 



and its solution yields 



2 arctan 



(11) 



where ipa is determined from the initial conditions. If the 
initial condition corresponds to the correlated pump and 
signal states (rj is non-zero), then ip = axcsm(jj/I) w 
rj/I. If signal and pump are initially not correlated, 
Po = and our approximation evidently fails because 
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FIG. 2: Schematic plots of the classically stable trajectories for the degenerate parametrical oscillator. 



all the time derivatives in Eqs. © stay zero all the time. 
In this important case we should go beyond the classical 
approximation and take into account quantum fluctua- 
tions of pseudospin components^. 

In order to do this we come back to the set of equations 
© but we will treat them as quantum equations for oper- 
ators (£, 77, Z). Using Eq. © we can see that at the initial 
moment, when I — \Z\ <C /, mean square fluctuations of 
£ and n are equal to (£ 2 ) = (rj 2 ) = N/4 + N/2{N/2- \Z\). 
At small £ and rj, \Z\ rj I and is approximately constant. 
We can therefore simplify Eqs. @ by putting Z = I, as 
follows 



V 



-Iri, 



(12) 



4V 

These equations are exactly solvable and give an expo- 
nential growth of r\ and exponential decrease of £. So, at 
longer time delays, when the approximate Eqs. Ip^l are 
not valid, we can still put £ « and consider the motion 
along the corresponding classically stable trajectory, but 
£, r] and Z should be treated as quantum-mechanical op- 
erators. The commutation relations at £ — * simplify: 



[V,Z) = 0, [S,Z]=ir), [ri,Z]=iZ. 



(13) 



Therefore, one can represent operators Z, r\ and £ as 
£ = —id /dip, Z — pcosip, r\ = pshup where p and <p are 
classical (not operatorial) variables. In this approxima- 
tion, the Hamiltonian (JSJ) corresponds to the following 
time-dependent Schroedinger equation: 



2Vp / . <9* 



4" cos if 



(14) 



This equation allows for explicit integration, namely 



* 2 sin^ = F (— -£-ln 



cp 

tan — 
2 



(15) 



where F is a function to be determined from the initial 
conditions. It is worth noting that |-F|<2y? = |W| 2 cZZ gives 
a probability distribution over Z for the degenerate para- 
metric oscillator. 

Eq. (|15fl shows that in the limit of large number of 
polaritons, N 3> 1, the dynamics of our system is es- 
sentially classical and its wavefunction can be replaced 
by the classical probability distribution, which, in turn, 
satisfies the Liouville equation^. Actually, comparing 
the arguments of the function F in Eq. (|15f) and clas- 
sical solutions for ip(t) given by Eq. Ijlll) . one can see 
that the probability distribution F(tp, t) at any time mo- 
ment t can be obtained from the initial one, ^0(^0)5 a s 
F(ip,t) = Fo((p(t)), where <p(i) is found by solving the 
classical system (0. In other words, one needs just to 
solve equations 10, treating Z, £, and r\ as classical vari- 
ables, for a set of different initial conditions with the 
statistics reflecting the initial \^\ 2 , and then the prob- 
abilities for all the pseudospin components at any time 
can be obtained. 

The validity of our approximation can be checked if we 
notice that on a short time-scale after an initial moment 
our system is equivalent to the harmonic oscillator with 



the inverted potential (Hint oc p 



where p and x 



are generalized momentum and coordinate respectively) . 
For such a problem the quasi-classical approximation is 
valid when the spread of the wavefunction is larger than 
the oscillator length (root mean square of x in parabolic 
potential)2i. In our case, when almost all the polari- 
tons are initially in the pump state, the wavefunction 
at t = is an eigenfunction of the operator Z. On the 
other hand, this function is the first eigenfunction of the 
harmonic oscillator. So, initially, the width of the wave- 
function is the same as the oscillator length, but in-plane 
components of the pseudopsin rapidly increase with time 
(see Eq.(12)), and the wavefunction spreads into the re- 
gion of validity of the quasi-classical model. This means 
that the dynamics of the pseudospin is quasi-classical at 
any time except the short initial time range, where the 
system obeys analytically solvable Eqs. (12). 

The above analysis refers to the first half-period of the 
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hyperspin motion. For the second half-period the same 
arguments can be used with the only difference that rj — > 
serves as the momentum operator and £ becomes the 
classical coordinate variable. 



B. General case 

The description of the non-degenerate parametric os- 
cillator (with different signal and idler states) is analo- 
gous to the previous case. Here, to deal with three po- 
lariton states, we introduce a 9-dimensional hyperspin 
pseudovector. Its components describe correlations be- 
tween all pairs of states. For example, "signal-pump" 
components are 

X ± = \ ( a s a t + a p a t) ( 16 ) 



Y x = -- (a s al-a p al) 



z i = 2 ( a l a p ~ a s fls ) 

The remaining components X 2 , Y 2 ,Z 2 describing corre- 
lations between pump and idler and X 3 ,Y 3 ,Z 3 for the 
idler-signal correlations can be obtained by respective 
changing of indices. As it can be readily seen from the 
definition, Z\ + Z 2 + Z 3 = 0. The commutation relations 
for the hyperspin components with the same index were 
already obtained in the previous section. Other compo- 
nents obey the following relations 



[X a , X a+ ±] — [Y a+ i, Y a ] — -Y a -i, (17) 



[X a , Y*+i] — [Y a ,X a+ i\ — —X a —i, 



[X a , Z a +i] — [X a ,Z a -i] — -^Yai 



[Y a , Z a +i] — [Y a , Z a —x] — — —X a , 



[Z a , Z a+ i] — 0. 

where a € {1, 2, 3}, a±3 = a. The following combination 
plays the role of the squared total angular momentum: 

9 i\t2 

J2(X? + Y? + ±Z*) = ?- + N. (18) 



The interaction Hamiltonian can be expressed in terms 
of the hyperspin operators in a very simple form 



Free propagation terms can be written as 
Wo = I [{E P - Ei)Z x + {Ei - E P )Z 2 + (E s - Ei)Z 3 - 



(20) 



{E s +E p + Ei 



N 



Our next step is to obtain a set of Heiscnberg equa- 
tions for the hyperspin operators. Using commutation 
relations 1|17H and simple algebra we come to 



W 

Z 1 = - — (X 1 Y 2 -Y 1 X 2 ), 
n 



(21) 



{X X Y 2 - YiX 2 ), 



V 



(XIY3 + Y1X3 + 2ZIY2 



En — E s 



V . 



Y3X2 — ZY1Z2 



X 2 = -j{X 3 Y 2 

Y x = —(YjYs - X X X 3 + 2Z X X 2 ) 
n 



Yo 



V 



(X 3 X 2 + Y 3 Y 2 + 2X X Z 2 ) 



V 



X 3 — — — (X1Y1 + Y1X1 — X2Y2 — Y2X2) 
n 



Yi, 
Y 2 , 
Xi, 

H X - 
Ei - E, 





h 




Ei 




E p 




h 




E p 




E s 




h 




Ei 




Ep 



Yo 



Y 3 = --(Xf-Y 1 2 



Ej, — En 



-x 3 . 



Further treatment of this system is very similar to one 
we applied in the degenerate case. The hyperspin system 
(|21|l has a classically stable trajectory defined by the fol- 
lowing conditions: 



Y% — X\ . X 2 



-Yi, Y X =±X X . 



(22) 



The same argumentation as in the Sec. IIA shows that 
in the limit of large number of polaritons the motion of 
system is concentrated near this trajectory. This fact 
allows us to reduce number of relevant variables. In order 
to do it, it is convenient to use another set of variables, 
namely 



H int = 2V(X 1 X 2 + Y 1 Y 2 ). 



(19) 



(X 1 +Y 2 ±X 2 TY 1 ), (23) 
(X 2 + Y 1 ±X lT Y2), 

C± = ^(z 1 -z 2 tx 3 ). 

On the classically stable trajectory defined by Eg. (|22fl . 
£,± = 0, rj± = X x =p Yj and C± = Z x =p X 3 /2. These vari- 
ables are remarkable by the fact that the commutation 
relations of the variables with the same index (+ or — ) 
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are the same as commutation rules for the operators of 
the angular momentum: 

[V± . £t] = iC± , [£±, C±] = "7± , [i?±, C±] = -iC± > 

(24) 

while the commutators of variables with different indices, 
and of any of them with Y3, are either exactly, or approxi- 
mately equal to zero in the vicinity of the classicaly stable 
trajectory Eq. I|22|) . When the hyperspin is close to the 
classically stable trajectory, the commutator of rj± and 
C± becomes negligible and the following substitution is 
possible 



(±=P±cosip±, rj± = p±sm(p±, 



(25) 



where p± and tp± can be treated as c- numbers. To satisfy 
first two relations in Eq. 124|) . one should take 



£± = -i 



d 



d<p± 



(26) 



taken into account phenomenologically within the hy- 
perspin approach, by adding linear dissipation terms to 
Eqs.J2U, m analogy to Bloch equations. In the most 
general case, to keep linear dependence of the decay rate 
on the population of polariton states, such terms should 
have the following form: the total polariton population 
(which has been so far assumed constant) is now given 
by the equation N = —aN — 2(3-1, where a and all 
the components of the vector f3 are constants; and the 
dissipation of the hyperspin vector is determined by the 
term — al + 51, where S is a second-rank tensor. 

(We remind that I is the 9-dimensional hyperspin pseu- 
dovector.) The parameters a, (3, and 5 should be chosen 
in such a way as to provide given values of the inverse 
lifetimes 7/. (k — s,p,i) of polaritons in signal (s), pump 
(p), and idler (i) states for any direction of the hyperspin. 
This condition is met if 



In terms of the new variables, the interaction Hamil- 
tonian Eq. (J2J) is decomposed in two parts describing in- 
dependent degenerate parametric oscillators 

,=± 7=± V ov i J 



■x = 1/3(7 S +7 P +7i), 
Pzi = 1/3(7 P - Is), 
Pz 3 = V 3 (7i -Hp), 



(29) 



(27) 

Here we assume that the parametric oscillator wavefunc- 
tion \& depends on p± and (p±. At the exact resonance, 
E p — E s = Ei — E p and free propagation term in the 
Hamiltonian may be omitted, since all possible configu- 
rations of the parametric oscillator have equal energies. 
The solution of the corresponding Schocdinger equation 
is a function of two arguments, namely 



"J/ 2 sin ip + sin <p_ 



(28) 



(3 Zs = 1/3(7 S - 7i), 
and the components of tensor 5 are 

Szi,z a = -5z 3 ,z 2 = Pz 2 , $z 1 ,z a = 



3Z2.Z3 



-Pz 3 



(30) 



fiz 2 ,Z! — —Sz 3 ,z-i — —fiz x , 
and all the other components of /3 and S are zero. 



'( 2 T t «- h 



tan ■ 



2Vp- 



t - In 



tan 



The function F is determined from initial conditions in 
the same way as for the two-level model. The physical 
meaning of Eq. (|28[) is exactly the same as in the case of 
the two-level model: we can see that, for a large number 
of polaritons, the wavefunction can be replaced by a clas- 
sical probability distribution which satisfies the Liouville 
theorem and evolves in accordance to the set of dynami- 
cal equations for hyperspin components (Eq.l|2iy). 



C. Finite polariton lifetimes 

The dynamical system considered above is idealysed 
as it does not include dissipation processes. In reality, 
polaritons have finite lifetimes dependent on the qual- 
ity factor of the cavity and on the non-radiative broad- 
ening of the exciton. The decay of polaritons can be 



III. NUMERICAL RESULTS AND DISCUSSION 

We will focus on the situation when the pump state 
was initially populated while signal and idler states are 
empty. First of all, we will consider the pump state with a 
definite number of particles (Fock state). The probability 
distribution for Xi and Yi components is given by 



(2^)% 2 a 3 



■ exp 



i*(*=o)r 



X\ + XI + Y? 



(31) 



2a, 



where 01 and 03 are the mean square values of the re- 
spective hyperspin components. The Z-projections of the 
hyperspin take their maximum values — ±iV/2, re- 
sulting in o- 3 = and o x = AT/ 4 (see Eq. JTSJl). 

Figure 3 shows the dynamics of the signal and pump 
populations calculated within our hyperspin formalism 
compared to those obtained by direct diagonalization of 
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FIG. 3: Signal (solid) and pump (dashed) populations as a 
function of time. Idler population coincides with signal one. 
Main figure corresponds to the results obtained in the hy- 
perspin formalism, and inset shows populations obtained by 
the direct diagonalization of the interaction hamiltonian. The 
difference in the behavior at time large than 7 ps may be at- 
tributed to the low (TV = fOO) number of particles in our 
simulations. 



if the initial value of <p+ is (po, then 
2Vp + 



In 


tan^ 


= - In 










tan T 



(32) 



(This relation can be derived comparing the argument of 
the function F in Eq. (|28|l at two time moments and 
t) . Thus the quarter of period at fixed (/3q is given by 



h/i = 



2Vp, 



In 



tan 



(33) 



Eq. (|33|l should be averaged with the initial distribution 
of i^o ■ Reminding that i^o < 1, one can put ipo ~ Vo/p+ 
where 770 is the initial value of rj + and p + = N/2. The 
variable rj+ = X± — Y\ has initially Gaussian distribution 
with the standard deviation <j v = 2a \ = N/2. Finally, 



T1/4 



x 



exp 



2a 7 



tl/4 = 



(34) 



the Hamiltonian Eq.© (see inset). The fully numerical 
solution of the Schoedinger equation with the Hamilto- 
nian given by Eq. (J2J) has been carried out by projecting 
the Hamiltonian operator to the basis of the Fock states 
with the definite number of particles in signal, pump and 
idler. The eigenvalues and eigenvectors of the Hamil- 
tonian matrix were found numerically making possible 
to construct a time-dependent solution in the form of a 
Fourier expansion. Using this time-dependent wavefunc- 
tion it is possible to calculate all relevant values, namely, 
the average populations of states and their second order 
coherence (see below). Within this approach, the total 
number of polaritons, N, cannot be taken very large. We 
performed calculations with N — 100, which provided 
reasonable computation time and, on the other hand, al- 
lowed the comparison with our quasiclassical hyperspin 
model valid at large N. Calculations within the pseu- 
dospin approach were performed by solving Eqs. (|21|l for 
a set of initial conditions distributed randomly according 
to Eq. (|31|l . with a subsequent averaging. We took inter- 
action constant V = 2.1 x 10 -2 meV so as to make the 
product VN ~ 2 meV close to the realistic valued. Pop- 
ulations of signal and pump states demonstrate damped 
oscillations centered at N/4 and N/2 respectively. Figure 
3 clearly demonstrates two main features of parametric 
oscillations: (a) the pump polaritons pass to signal and 
idler states and vice versa, and (b) polaritons arrive to 
the signal state with some delay with respect to the ini- 
tial moment when the pump state was excited, as it was 
shown experimentally in Ref. lift 

The quarter-period of the oscillations T1/4 can be 
found analytically from the solution of the Schroedinger 
equation obtained in a previous section. For simplicity 
we consider an evolution of variables Lp + and £ + only, as- 
suming that the initial distributions in + and — variable 
sets are the same. According to Eq. (|2*5|) . the quarter of 
the period corresponds to (p + = n/2. On the other hand, 



h 

VN 



In 



N 



0.9 ps 



where 7 « 0.577216 is the Euler constant. This value 
concides with the result of numerical calculation ^1/4 ss 
0.9 ps (see Fig. 3). 

The damping of the oscillations results from the fact 
that the Hamiltonian @ has almost continuous spec- 
trum at N ^> 1. The summation of oscillations on 
different eigen-frequencies leads to the degradation and 
damping of the oscillations of the signal, idler and pump 
populations. The steady state values can be obtained 
from the detailed balance, i.e. the numbers of incoming 
and outcoming particles for each state should be equal. 
Since polaritons can arrive to or leave the pump state 
by pairs only (see Fig. la) and both signal and idler 
were not populated initially, then the steady state values 
for the pump, N p , and signal, N s , states obey the equa- 
tions N p + 2N S = N and 2N S = N p with the solutions 
N p = N/2 and N s = N/4. 

One can see that at time t < 7 ps the hyperspin for- 
malism and direct diagonalization give utterly similar re- 
sults. The discrepancies at longer time delays are due to 
the fact that the number of particles (N — 100) we used 
for simulation was too small for the quasi-classical ap- 
proximation to hold on the larger time scale. 

Now we proceed to discussion of quantum statistics of 
the parametric oscillator. We consider two different cases 
of initial pumping: coherent with initial pump statis- 
tics P(N P ) = a 2Np exp (— a 2 )/N p \, where a 2 is the av- 
erage number of particles and thermal one with initial 
statistics P(N p ) = (1 - 9)6 N " , where 9 = N p /(1 + N p ) 
and N p is the average number of polaritons in the pump 
state. Figure 4 presents time dependence of the signal 
and pump populations and so-called second-orded coher- 
ence g 2 = (aj- 2 afe)/(oJ.afc) 2 (with k = s,p) which can be 
measured in the two-photons counting experiments. Its 
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Time (ps) Time (ps) 



FIG. 4: Time dependence of the signal (solid) and pump (dashed) populations in the case of coherent (a) and thermal (b) 
pumping without lifetime. Inset shows second-order coherence for the signal state by solid line and for the pump state by 
dashed line. 




Time (ps) Time (ps) 



FIG. 5: Time dependence of the signal (solid) and pump (dashed) populations in the case of coherent (a) and thermal (b) 
pumping with lifetime 5 ps. Inset shows second-order coherence for the signal state by solid line and for the pump state by 
dashed line. 



value characterizes the statistics of the relevant state and 
equals to 1 for the coherent state and to 2 for the thermal 
one. It can be seen that in the case of coherent initial 
pump the time dependences of the populations are simi- 
lar to those obtained for the pump in the Fock state (see 
Fig. 4(a) and Fig. 3 for comparison). The case of ther- 
mal pumping (Fig. 4(b)) is different. The oscillations 
are suppressed due to the large spread of the initial dis- 
tribution in the particle number and the populations of 
the states reach their steady values almost monotonously. 
Our calculations show that the steady state values of gl 
and <7p are very close in magnitude and are 1.7 in the 
case of initially coherent pump and 3.3 for the thermal 
pump. These values are confirmed by the direct diago- 
nalization of the hamiltonian (N — 100) which gives 1.7 
and 3.2 respectively. Such values are different from those 
obtained in Ref. [15| for constant wave pumping. It may 
seem strange that coherence dissapears in the interacting 
system. However this result becomes clear if one notices 
that in the case of the polariton laser the spontaneous 



coherence build-up is possible if the incoming scatter- 
ing rate to the ground state is larger than the outgoing 
scattering rate^S. In our case of a three level system, in- 
coming and outgoing rates are actually the same and the 
coherence is destroyed by multiple scattering processes. 

The calculated populations and statistics for the po- 
laritons with finite lifetime arc presented in Figure 5. The 
lifetimes for all states were choosen identical and equal 
to 7 _1 = 5 ps. One can see in the case of initially coher- 
ent pump (Fig. 5a) that the populations of signal and 
pump states behave nonmonotonously. From the very 
beginning polaritons start passing from the pump state 
to the signal and idler states, while initially this process 
goes spontaneously, and signal population growth slowly. 
Then, stimulated scattering switches on, and both signal 
and pump populations change rapidly. Then the process 
of return of polaritons from signal and idler to the pump 
state starts, etc. Because of the finite lifetime, the total 
number of particles in the system decreases. When it falls 
below the stimulation threshold N ~ H/Vj no more os- 
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cillations of the occupation numbers of signal and pump 
states can be seen. 

The situation is different in the case of the thermal 
pump (see Fig. 5b): no pronounced oscillations can be 
seen for the pump population, however the sharp bend 
of the time dependence of the pump population and the 
increase of the signal population demonstrate that the 
parametric process takes place nevertheless. 

Interestingly, introduction of lifetime in the system re- 
sults in better coherence at t — 14 ps g% w 1 in the case 
of coherent pumping and in almost the same value g% ss 3 
for the initially thermal pump. 

In conclusion, we have presented the general formal- 
ism describing the dynamics of the optical parametric 
oscillator based on a semiconductor microcavity in the 
strong coupling regime in the case of large number of po- 
laritons. The introduction of the hyperspin allows to ob- 
tain a quasi-classical solution of the quantum parametric- 
oscillator problem. We have shown, that the probability 
distribution for the hyperspin components obeys the Li- 
ouville equation. Our approach is shown to give very 
good agreement with the method based on the direct di- 
agonalization of the Hamiltonian for a "moderate" num- 
ber of particles (N = 100) allowing such a numerical 
solution within a reasonable computation time. At the 



same time, our quasiclassical hypespin model radically 
reduces the computational complexity of the problem for 
large number of particles, gives a possibility to introduce 
finite lifetimes of polaritons, and allows to obtain some 
analytical results such as an expression for the period 
of oscillations. The presented numerical results suggest 
that the polariton parametric oscillator can be used as a 
system to grow coherence in the case where irreversible 
processes are present. Our formalism can be extended 
to allow for spin degree of freedom of polaritons and 
might therefore be suitable for description of polarization 
properties of the microcavity emission in the parametric 
regime. 
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